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ABSTRACT 

We examine 14 hot white dwarfs for signs of infrared excess using the Spitzer 
Space Telescope. Twelve of the objects are massive white dwarfs which have 
been suggested to be the result of binary mergers. The remaining two objects 
are undermassive white dwarfs which again may be the result of mergers or the 
inspiral of a substellar companion. In no case do we find any evidence for signifi- 
cant infrared excesses out to wavelengths of 8/im. This places strong constraints 
on the presence of orbiting dust and weaker constraints on the presence of close 
substellar companions. 

Subject headings: white dwarfs - planetary systems: protoplanetary disks - in- 
frared: stars 



1. Introduction 

The mass distribution of white dwarf stars is one of the fundamental measures of stellar 
evolution. The majority of white dwarfs have masses distributed in a narrow range between 
0.5 - 0.6M© (Weidemann & Koester 1984; Bergeron, Saffer & Liebert 1992 -BSL; Liebert, 
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Bergeron & Holberg 2005 -LBH). As studies of this type grew steadily more accurate, two 
additional, smaller peaks appeared in the distribution (BSL, LBH). The first is centered 
around a mass ~ O.4M and represents the outcome for stars in close binaries, whose stellar 
evolution has been truncated by Roche lobe overflow before they reached the Helium flash. 
The second minor peak lies at higher masses > 0.8M Q . This excess of massive white dwarfs is 
even more extreme amongst white dwarfs selected from ultraviolet or X-ray surveys (Marsh 
et al 1997; Vennes et al 1997) because surveys biased towards hotter stars are more likely to 
find massive white dwarfs (a result of their smaller radii). 

This apparent excess of massive white dwarfs is of interest because it represents a poten- 
tial sample of white dwarf merger remnants. A natural outcome of stellar binary evolution 
is a significant population of double degenerate white dwarf binaries in close orbits, close 
enough that they should merge within a Hubble time (Webbink 1984; Iben 1990; Yungelson 
et al 1994). Most of these binaries have total masses well below the Chandrasekhar mass 
limit and so are expected to form a single isolated remnant star. The expected population of 
> O.8M degenerate merger remnants is thus quite naturally associated with the observed 
excess of massive white dwarfs. However, there is as yet no convincing evidence that any of 
these massive white dwarfs are indeed the result of a merger. 

Numerical simulations of the merger of two unequal mass white dwarfs suggest that the 
less massive of the pair is tidally disrupted and accreted onto the more massive one in the 
form of a thick accretion disk (Mochkovitch & Livio 1989; Benz et al 1990; Guerrero, Garcia- 
Berro & Isern 2004). Thus, it would appear that the natural outcome of such a picture is 
a hot white dwarf surrounded by a thin remnant disk which contains most of the angular 
momentum. The nature of this disk is likely to be very unusual, given that white dwarfs are 
composed primarily of Carbon and Oxygen, rather than the Hydrogen and Helium of more 
traditional circumstellar disks. The observational signature of such circumstellar disks is an 
infra-red excess, as dust particles in the disk absorb stellar photons and re-radiate them at 
longer wavelengths characteristic of the cooler temperatures in the disk. Indeed, just such 
an excess is observed for the white dwarf G29-38 (Zuckerman & Becklin 1987), although 
the source of the reprocessing disk is likely to be something else in that case (Graham et al 
1990; Debes & Sigurdsson 2002; Jura 2003). Thus, as a test of the merger hypothesis for the 
origin of massive white dwarfs, we have undertaken a Spitzer program to search for infra-red 
excesses around the youngest /hottest massive white dwarfs. 
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2. Observations 

The sample was selected from those white dwarfs found in the EUVE and ROSAT 
surveys, which were optically classified as being massive (> 1M ). Selection from these 
surveys guarantees that the white dwarfs are hot and therefore young (effective temperatures 
> 35, 000K correspond to ages < 10 7 years). This strategy has its strengths and weaknesses. 
It is observed that infra-red excesses around main sequence stars disappear on timescales 
of a few xlO 8 years (Habing et al 1999; Rieke et al 2005, and references therein), so that 
we placed a great emphasis on observing stars in the earliest possible stages of their present 
incarnation. On the other hand, LBH suggest that, while there is definitely an excess in 
massive white dwarfs over that expected from single star evolution, the UV/X-ray samples 
may be biased towards the progeny of nearby O/B associations which may dilute the fraction 
of merger remnants. We shall return to this subtle point in § 5. 

Drawing from the EUVE (Vennes et al 1997) and ROSAT (Marsh et al 1997) samples, we 
chose twelve white dwarfs with masses > 1M (based on the optical spectral classifications). 
To this sample of potential merger remnants we also added two low-mass white dwarfs which 
have also been suggested as merger remnants (Maxted & Marsh 1998) because they show 
no evidence for a companion. Table 1 shows the resulting sample of fourteen objects, with 
temperatures, masses and distances. 

In order to determine the appropriate bandpass, let us consider briefly the kind of 
remnant disk one might expect from a merger remnant. White dwarfs are more compact 
than main sequence stars and the angular momentum deposited in the disk is only sufficient 
to produce disks of extent ~ 10 11 cm. As such, we expect the disks to be somewhat hotter than 
traditional protoplanetary disks. Simple models suggest that the emission from a passively 
reprocessing disk will peak at wavelengths ~ 10/xm. Thus, observations were planned in the 
Spitzer IRAC bandpasses. The spectral coverage of Spitzer from 3.6-8/xm is well-suited to 
the expected spectral shape. These issues will be discussed further in § 4. 

The Spitzer observations were taken in all four IRAC bandpasses, with 3 frames of 12s 
exposures in a 5 dither pattern. The resultant total exposures were 180s for each band. The 
exposure lengths were chosen to produce detections of simple reprocessing models at levels 
~ lOOo". The exposures were also sufficiently long to ensure confident detection of the white 
dwarf photospheric emission at 3.6/xm, allowing us to place stringent constraints in the event 
of no detection. 
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3. Results 

Table 2 shows the observed fluxes. We include V-band magnitudes from the literature 
and have either measured 2MASS J-band fluxes or culled results from the literature. There 
is no convincing evidence of an infra-red excess in any of the sources. In several cases, the 
observations are good enough to trace the photospheric emission all the way to 8/im and in 
all cases the photosphere is detected at expected levels in the shortest IRAC bandpass. 

Matching finding charts between optical and IRAC bands can prove challenging, espe- 
cially for objects as blue as ours. Thus, we provide IRAC finding charts for our objects in 
Figures 1-7. 

3.1. Noteworthy Objects 

3.1.1. GD 50 

The prototypical supermassive white dwarf, GD 50 was classified as such by Bergeron et 
al (1991) on the basis of the anomalous width of the Balmer lines in the spectrum. Bergeron 
et al derived T e ff = 43500 ±1500 K and log g = 9.00 ±0.15 (both 3<x errors), yielding a mass 
~ 1.2M . They noted that, under the assumption that massive white dwarfs are born from 
massive stars, one would expect the total stellar age of this star to be short (~ 10 7 years) 
so that it was rather surprising to find no evidence of recent star formation in the vicinity. 
Following this line of argument, one may consider GD 50 one of the most hopeful candidates 
for finding evidence of a recent merger. 

The white dwarf is solidly detected in the first three IRAC bandpasses and marginally 
detected in the fourth bandpass. The resulting photometry is shown in Figure 8. We 
see that although the object is detected in the optical, in 2MASS and by Spitzer, all flux 
measurements are consistent with photospheric emission. We show a 40,000 K black body 
as a solid line. Matching the model fluxes to the observed fluxes, assuming the spectroscopic 
temperature, provides a constraint on the ratio of radius (Rg in units of 10 9 cm) to distance 
(D pc in units of parsecs) 

R 

— - = 0.0108 ±0.0016 

which is in excellent agreement with the values quoted in Bergeron et al (R = O.OO57i? and 
D ~ 37pc, which imply Rg/D pc = 0.011). Thus, our measurements of GD 50 are consistent 
with its status as a massive white dwarf if the distance is of order the nominal 37 parsecs. 
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3.1.2. EUVE J0916-19.7 



Green et al (2000) claim a significant near-infrared excess for this object, with J and K 
values approximately 4 magnitudes above the values predicted from atmosphere models. As 
can be seen from Figure 9, they have almost certainly reported the magnitude for the star 
~ 7.2" away from the true source. The circled white dwarf has flux measurements consistent 
with photospheric values. To confirm this, Figure 10 shows an I-band image of the field taken 
with the Palomar 60" telescope. The white dwarf is circled. Figure 11 shows the near-IR 
fluxes for both the white dwarf and the nearby star. 

We can determine R 9 /D pc for these objects as well. The red object is probably a main 
sequence star at some distance behind the white dwarf. 



This object has been suggested as a merger remnant on the basis of several lines of 
argument. Not only is it massive, but it is also magnetised and has a rapid rotation period 
of 725s. All these features can be plausibly explained by a double degenerate merger (Barstow 
et al 1995; Ferrario et al 1997). 

Figure 12 shows the resulting field. Once again the object is easily detected in the first 
three bands and marginally detected in the longest bandpass. However, all the detections 
are consistent with photospheric emission from a 40,000 K white dwarf. The results are 
shown in Figure 13 along with the flux measurements for the close white dwarf companion 
LB9802. This white dwarf has a spectroscopic temperature and gravity (T e ff = 16030 ±230, 
log g = 8.19 ± 0.05 - Barstow et al 1995). It was reported as being at a separation of 16" by 
Barstow et al, but the Spitzer image yields a separation of 7.4" (a similar number is obtained 
from 2MASS). 

The temperature determination for EUVE J0317-85.3 is T eff ~ 40,000 K, but there is 
no spectroscopic gravity. However, Barstow et al argued for the binarity of the two white 
dwarfs (on the basis that their close separation was unlikely to be random). Placing the 
two objects at the same distance then implies a high mass for the hotter star because it is 
fainter. Our Spitzer photometry yields a similar answer. In the Rayleigh- Jeans tail, the flux 
ratio (assuming the same distance) is 



if we designate LB 9802 as star 1 and EUVE J0317-85.3 as star 2. Using the measured 



3.1.3. REJ03 17-853 
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temperatures, this implies R1/R2 ~ 2.5. So, for R\ ~ 8 x 10 8 cm (based on the measured 
gravity of LB9802), this yields R 2 ~ 3.2 x 10 8 cm. 

3.1.4. WD0235-12.5 

Figure 14 shows the field around WD 0235-125. Although once again no credible excess 
is detected (see Figure 15), there is an additional source within 4" of the white dwarf. Any 
claim about close companions in fields as crowded as these must be viewed with caution, 
but the characteristics of this object are of interest. The fluxes shown in Figure 15 are 
consistent with a blackbody of ~ 5000K and radius approximately twice as large as that of 
WD 0235-125 (if the two are physically associated). The separation of 3.6" at a systemic 
distance of 6.6 pc means an orbital separation of 238 Au. This would make it the second 
star in our sample to have a white dwarf companion. However, the cooling age of such a 
companion is ~ 3 x 10 9 years (as opposed to ~ 10 8 years for WD 0235-125). Perhaps the 
massive white dwarf is the product of mass transfer and binary evolution. One problem with 
this speculation is that a main sequence star with lifetime ~ 3 Gyr has a mass ~ 1.4M Q , so 
that the final white dwarf is of comparable mass to the putative progenitor. 

3.1.5. WD1614+136 & WD1353+409 

These two objects are different from the others in our sample. Identified as low mass 
Helium core white dwarfs, these stars are expected to have a close, dark companion which 
is responsible for the stripping of the progenitor envelope before the core mass got large 
enough to undergo Helium burning. Yet Marsh, Dhillon & Duck (1995) found no evidence 
for the radial velocity variations that characterise most of the stars in this class. Possible 
alternatives to the traditional route involve either the merger of two low mass white dwarfs 
(although Marsh & Maxted 1998 find normal rotation velocities) or perhaps stripping by a 
planetary mass companion (Nelemans & Tauris 1998). 

Figure 16 shows the fluxes measured for WD1614+136 and WD1353+409. Once again 
the emission is completely consistent with the photospheric expectations, from optical wave- 
lengths to 8/im. The field of WD1353+409 (Figure 17) shows a marginal detection of a 
nearby source that may be a companion. If it is indeed physically associated, then the flux 
ratio is consistent with a 5000 K white dwarf companion. Larger radius objects such as brown 
dwarfs or planets require lower temperatures which are inconsistent with the photometry (see 
§4.3.1). 
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4. 



Models 



Given that all of our observations are consistent with photospheric emission, what upper 
limits can we place on the presence of a disk? 

Despite the fact that conditions in merger remnant disks are far from the usual situation, 
some models do exist (Livio, Pringle & Saffer 1992; Menou, Perna & Hernquist 2001; Livio, 
Pringle & Wood 2005). For metallic materials, once the disk temperatures decrease below 
~ 10 4 K the material will recombine and one might expect the formation of dust. Simple 
models of viscous disk evolution in a merger scenario (Menou et al 2001) suggest that the 
outer edge of the disk can evolve to scales ~ 4 x 10 n cm before the entire disk becomes 
neutral. Furthermore, setting the temperature of the passive disk equal to 10 4 K yields the 
inner edge of any plausible disk (interior to this the stellar irradiation can keep the material 
ionized even in the absence of viscous heating). This yields a value 



The simplest disk is a flat, optically thick, passively reprocessing disk in which each 
annulus radiates as a black body (Adams et al 1987 and references therein). We note that 
this is a more conservative model than that of Livio et al (2005), which flares as a result of 
vertical hydrostatic equilibrium. So our constraints on that model would be even stronger 
than the ones derived below. 

In the flat disk model, the disk luminosity is 



where a« and a Q are the inner and outer radii of the disk and i is the inclination angle of 
the line of sight. The temperature profile T(a) is determined by equating the incident and 
radiated luminosity at each annulus, 




Thus a plausible range of disk scales ranges from ~4x 10 9 



'cm to 4 x 10 cm. 



4.1. Flat Black Body Disk 
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where R* and T* are the stellar radius and effective temperature respectively. This results 
in an expected flux 

^ = 17mJy (^) (iJbr) (lo^) cosi i i^T 

for a disk around a star at 100 pc distance (z = hc/\kT). We have chosen numbers appro- 
priate to hot white dwarfs rather than the usual main sequence star parameters. We note 
that the only tunable parameters for this simple model are the inner and outer radius and 
the inclination. 

As can be seen from Table 2, our observations are all at the level of a fraction of a mJy 
so well below the predicted level. Essentially all of the reasonable parameter space is ruled 
out. Figure 18 shows an example of the model comparisons with the data for GD 50. In 
order to fit the data at the faint end, the inner edge of the disk has to be placed at ~ 10 11 cm 
(close to the outer edge expected in the models of Menou et al) and even then only fits for 
disks seen almost edge-on 1 . Thus, passive black body models are ruled out. Table 3 shows 
the constraint on cosi we infer by requiring that the model flux be less than the detected 
8/xm level (assuming Oj = 10 11 cm). 



4.2. Dust only disks 

Chiang & Goldreich (1997) introduced an extension to the standard disk model incor- 
porating a superheated dust layer in addition to the gaseous disk. The temperature balance 
for the dust is different because the emissivities in the UV/optical (where flux is absorbed) 
and in the infrared (where flux is emitted) are markedly different for dust. Thus, the disk 
can behave as if it is optically thick in absorption but optically thin in emission. 

Given the fact that the disk in the situation under consideration can be almost en- 
tirely metallic, it is of interest to consider the effects of a flat, dust-only, reprocessing disk. 
Modelling this disk in the manner of Chiang & Goldreich, we infer a different temperature 
distribution from the previous section 

where we have modelled the ratio of infrared to optical emissivities as oc (T/T*)' 3 . We will 
assume (3 ~ 1 as in Chiang & Goldreich. The emergent emission must also be diluted by a 



1 Furthermore, even the edge-on case may be ruled out in the case of moderate flaring as the disk would 
be expected to eclipse the star and lower the optical luminosity - Livio, Pringle & Wood 2005. 
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factor 1 — exp(—Tm) where tjr is the optical depth of the disk in the infrared. The resulting 
temperatures are somewhat higher at a given radius than in the blackbody disk and this 
model is shown in Figure 19, as the short dashed line. In this instance, varying the inner 
disk edge does not change the flux profile as much as before, so we show only the case of 
cij = 10 10 cm. This model is also ruled out as it requires all the observed systems to be 
essentially edge on in order to meet the constraints. 

If one reduces the mass of dust in the disk still further, the disk becomes optically thin 
even in the ultra-violet. In this case, the absorbed flux is reduced by a factor 1 — exp(— Tuv), 
where Tuv is the average optical depth of the disk to incident stellar photons. This additional 
parameter allows us to fix the disk inclination to more reasonable values and determine the 
level of t uv needed to fit the data. The long dashed curve in Figure 19 shows this optically 
thin disk with the same inner and outer edges as before, but now with an inclination of 45° 
and Tuv = 0.05. 

We can infer an approximate surface density £ ~ tuv/^uv where kuv is the appropriate 
opacity. For k uv ~ 10 3 cm 2 g~ 1 , this means £ ~ 5 x 10~ 5 g cm~ 2 , which yields an approximate 
disk mass 

M disk ~ nR 2 Z ~ 1.6 x I0 18 g < 10~ 15 M Q . 

This is the mass of an ordinary asteroid in the solar system - a very stringent limit on the 
amount of dust remaining in the system. 

Table 3 also shows some quantitative limits for each system in the optically thin limit. 
We place constraints on tuv assuming an optically thin, superheated dust disk inclined at 
60°. One can also translate that into an estimate of the dust mass. 

4.3. Planets 

We have thus far placed constraints on passively reprocessing dust disks. If the material 
forms rapidly into planets, then the appropriate signature is not that of a dust disk, but 
of a young, hot planet. However, such planets will be almost impossible to detect. Since 
such planets would be made of the same material as a white dwarf they obey roughly the 
same cooling physics. Indeed, the young ages of the systems under discussion mean that 
the physics is well approximated by the simple Mestel model (Mestel 1952; see D'Antona & 
Mazzitelli 1990 or Hansen 2004 for more extensive reviews). In that model, luminosity at 
fixed age is proportional to mass, so that the low mass planets will fade rapidly (much more 
rapidly than a gaseous planet of similar mass) compared to the white dwarf and will never 
become visible as an excess. 
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However, in the case of our two low-mass, Helium-core white dwarfs, there is a possibility 
of finding brown dwarfs or normal planets as companions. 

4-3.1. Low Mass White Dwarfs 

Two of our objects, WD1614+136 and WD1353+409, are low mass white dwarfs without 
any obvious companion. It has been suggested that they are still the product of truncated 
stellar evolution, but that the Roche Lobe overflow was induced by a planet or a brown 
dwarf (Nelemans & Tauris 1998). In this case, we might expect the signature of a more 
traditional planetary mass object. Figure 20 shows the effect on the photometry of a 10Mj 
brown dwarf of age 10 8 years (model taken from Burrows, Sudarsky & Lunine 2003 - BSL). 
This model is ruled out by the measurement in the IRAC2 band, which lies several a below 
the predicted level. Lower mass or older objects are consistent with this measurement. 

Figure 21 shows the various BSL models in terms of age and mass. Solid points indicate 
models ruled out as companions to WD1614+136 (in all cases it is because of an excess of 
4.5/j.m flux). Open circles are models which are consistent with the observations. We can 
rule out models down to 5Mj for very young ages (10 8 years) but cannot rule out models as 
massive as 25Mj for ages > 10 9 years. 

Similar constraints are found for unresolved companions to WD 1353+409. If we consider 
the nearby (~ 5") object to be a physical companion, then the flux ratio in the Rayleigh- 
Jeans tail is ~ 4 at short IRAC wavelengths. This is consistent with the neighbour being 
a white dwarf of similar radius but lower temperature (~ 5000-6000 K), in a ~ 630 Au 
orbit. Increasing the radius to that of a brown dwarf would imply temperatures ~ 100 K, 
which is not consistent with the observed photometry (for either a black body or one of 
the BSL models). Of course, it could always be a background G or K dwarf of appropriate 
temperature. However, for this line of sight, it would have to lie ~ 7 kpc above the Galactic 
plane. 

5. Discussion 

We have searched twelve hot, massive white dwarfs and two hot, undermassive ones, 
for signs of infrared excesses indicative of a circumstellar dust disk reprocessing stellar ra- 
diation. We have found no convincing evidence of any such excess. In all cases, the objects 
are detected to wavelengths as long as A = 4.5/im and are found to be consistent with pho- 
tospheric emission as predicted on the basis of their optically determined temperatures and 
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luminosities. In several cases this conclusion can be extended out to 8pm. 

Our negative result can be used to place quantitative constraints on disk models. For an 
optically thick, passively reprocessing blackbody disk, the observations place stringent limits 
- requiring either than the inclinations be nearly edge on in all cases or that the inner radius 
of the disk be located at distance > 10 11 cm (which is unlikely within the angular momentum 
constraints of the formation scenario). If we postulate a gas- free dust disk which is optically 
thin, we can place constraints on the amount of dust mass present. In all cases the mass 
cannot be larger than that of a solar system asteroid - and most require M < 10 19 g. 

Our motivation for this search was that the mass distribution of white dwarfs shows an 
excess of systems with masses > 1M Q . It has been suggested that such an excess could result 
from the mergers of two lower-mass white dwarfs driven by gravitational wave emission. If 
such a merger leaves behind a circumstellar disk of tidally disrupted material, one potential 
signature would be an infra-red excess resulting from reprocessed stellar radiation. 

Our stringent limits suggest that, if such disks do indeed form, the dust must be re- 
moved or evolve on timescales shorter than the cooling times of our sample white dwarfs 
(~ 10 7 years). There are several possible mechanisms. 

A dust disk does not remain in orbit indefinitely. Small grains are blown out of the sys- 
tem by radiation pressure and larger grains spiral into the star under the action of Poynting- 
Robertson drag. For our white dwarfs, radiation pressure removes only the smallest grains, 
because even for our hottest objects, L ~ 0.1L & . Poynting- Robertson drag is significant, 
especially because our disks are somewhat more compact than traditional protoplanetary 
disks. The inspiral time for dust under the action of this drag is 

_ 4tt < s > p s c 2 a 2 

1 p» — 

R 3 U 

where < s > and p s are the average grain size and density respectively, and a is the distance 
from the star, of luminosity L*. Numerically, this is 

T PR = 100 years (^) {^^) (wk^) " 

So, the removal time of dust in this system is very short. Nevertheless, a detectable signal 
is still expected if the dust is continuously replenished by asteroid or comet collisions in the 
debris disk. Simple estimates of collisional evolution in disks (e.g. Dominik & Decin 2003) 
suggest that collisions between asteroids occur on a characteristic timescale 

T cM = 40 years [^-) 



- 12 - 



for a disk of total mass M disk composed of asteroids of size R. Thus, our expectations 
correspond to a scaled-down version of a more traditional protoplanetary debris system 
surrounding a main sequence star. 

Nevertheless, we see no evidence of a dust signature - if there was a merger the dust was 
either removed very quickly before it could form larger bodies, or else it formed sufficiently 
few large bodies that the collision time is long and the dust replenishment is low. Requiring 
therefore that T co u > Tp^ and setting R = 4km (so that the limits on the observable excess 
correspond to the mass of a single disrupted asteroid), places a limit on the characteristic 
total disk mass in this simple model M^sk < 10 21 g. One can potentially raise this limit 
by a factor ~ 1/Z ~ 50 if one postulates that only refractory dust can survive at the 
temperatures inferred for these compact disks. If most of the Carbon and Oxygen that 
emerges from the disrupted white dwarf forms ices which are rapidly sublimated, then the 
amount of dust is limited by the available refractory elements. However, the resulting limit 
is still no larger than an asteroid. 

The fact that disks can indeed occur around white dwarfs and that they can be detected 
is shown by the infra-red excesses of the stars G 29-38 (Zuckerman & Becklin 1987) and 
GD 362 (Becklin et al 2005; Kilic et al 2005). These white dwarfs are also DAZ stars, 
which show evidence for trace amounts of Calcium and other refractory elements in their 
atmospheres (Zuckerman & Becklin 1987; Gianninis et al 2004). The disks in these cases 
are thought to be the products of the disruption of an asteroid (e.g. Jura 2003) rather than 
debris from a merger (while GD 362 is indeed quite massive, G29-38 has a very ordinary 
mass) . Nevertheless, this demonstrates that the basic model we set out to test is viable and 
that our null result is significant. 

It is, of course, also quite possible that there are one or more flaws in the merger model 
for the origin of massive white dwarfs. One possibility is if one of the two merging stars has 
a helium core. If this were the case, then the accretion of the helium onto the surface of the 
merger product produces a helium burning star rather than a star+disk system (Iben 1990), 
However, this eventuality seems unlikely for the massive systems considered here, because 
that requires that one of the components be < O.4M and therefore the primary would have 
to be massive to begin with. 

Finally, the massive white dwarfs may simply be the result of single star evolution after 
all. LBH confirm a significant excess of massive white dwarfs in the sample of older, cooler 
white dwarfs drawn from the Palomar Green (PG) sample (Green, Schmidt & Liebert 1986). 
Indeed, they assert that 80% of the white dwarfs in the PG sample with mass > O.8M 
are formed by some other process other than single star evolution. Nevertheless, they also 
note that the fraction of massive (> O.8M ) white dwarfs with M > 1M is higher in the 
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UV/X-ray selected surveys than in the PG survey. This might be due to the fact that the 
former are selected from an all-sky survey while the latter is restricted to high Galactic 
latitudes. LBH suggest that the EUVE samples may have a higher fraction of massive white 
dwarfs born from single star evolution than the PG sample because the O/B stars in Gould's 
belt may provide a significant enhancement in the current local birthrate of such objects. 
Thus, it is still not entirely clear how much of a binary contribution is expected to contribute 
to our sample (drawn from the EUV/X-Ray catalogues to get the hottest and so youngest 
systems). Finally, the constraints on the true shape of the initial-final mass relation for 
stars are still rather weak and it is possible (although far from proven) that the features in 
the mass distribution are simply the result of a non-linear relation (e.g. Ferrario et al 2005) 
between the main sequence and white dwarf masses of stars. 

In conclusion, we have examined 14 hot white dwarfs for signs of near-infrared excess. 
These have all been suggested as potential merger remnants (12 supermassive white dwarfs 
and 2 undermassive white dwarfs). In no case do we find any evidence of infrared excess. 

This work is based on observations made with the Spitzer Space Telescope, which is 
operated by the Jet Propulsion Laboratory (JPL), California Institute of Technology (CIT) 
under National Aeronautics and Space Administration (NASA) contract 1407. We thank 
NASA, JPL and the Spitzer Science Center for support through Spitzer contracts 1264152 to 
UCLA. BH acknowledges the support of an Alfred P. Sloan Foundation Fellowship and thanks 
Mike Jura, Jim Liebert, Ben Zuckerman and Jay Farihi for comments. This publication 
makes use of data products from the Two-Micron All Sky Survey, which is a joint project 
of the University of Massachusetts and the Infrared Processing and Analysis Center/CIT, 
funded by NASA and the National Science Foundation. This publication has also made use 
of the online McCook & Sion White Dwarf Catalog and NASA's Astrophysics Data System. 
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file. Figure fl.jpg goes here 

Fig. I. — The circle indicates the location of the white dwarf in each of these fields. Both 
images were taken in the IRAC 2 band. Each image is approximately 6' x 6.4'. 

Figure f2.jpg goes here 

Fig. 2. — The circle indicates the location of the white dwarf in each of these fields. 

Figure f3.jpg goes here 
Fig. 3. — The circle indicates the location of the white dwarf in each of these fields. 

Figure f4.jpg goes here 
Fig. 4. — The circle indicates the location of the white dwarf in each of the fields. 

figure f5.jpg goes here 
Fig. 5. — The circle indicates the location of the white dwarf in each of the fields. 

figure f6.jpg goes here 
Fig. 6. — The circle indicates the location of the white dwarf in each of the fields. 

figure f7.jpg goes here 



Fig. 7. — The circle indicates the location of the white dwarf in each of the fields. 
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GD 50 




2 4 6 8 



Fig. 8. — The points show flux measurements in the four IRAC bands, the three 2MASS 
bands and in the V band. The solid line is for a 40,000 K black body. There is no evidence 
of an infrared excess. 

Figure f9.jpg goes here 

Fig. 9. — The circle indicates the location of the white dwarf EUVE0916-19.7 in the IRAC 1 
bandpass. The star to the upper left is most likely the source of the infra-red excess reported 
by Green et al., as shown in Figure 11 



Figure flO.jpg goes here 



Fig. 10. — The circle indicates the location of the white dwarf EUVE0916-19.7. This is an 
I-band image taken with the Palomar 60" telescope. 




Fig. 11. — The lower set of points are the flux measurements (4 IRAC bands and the BVRI 
magnitudes) for EUVE J0917-19.7 (the circled source in Figure 9). The two solid circles near 
the top left are the J and K fluxes Green et al claim to measure for EUVE J0917-19.7. The 
four open circles are the IRAC fluxes of the object 7.2" away from the center of the circle in 
Figure 9. The dashed curve is a black body of temperature 8000 K, while the solid curve is 
for a blackbody of 56,400 K. 
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Figure fl2.jpg goes here 

Fig. 12.— The circle shows the location of the white dwarf REJ0317-853 in the IRAC 2 
passband. The star to the upper left is also a white dwarf, LB 9802. 



EUVE J0317-85.5 
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Fig. 13. — The solid points and curve are the flux measurements for EUVE J0317-85.3 while 
the open points and dotted curve are for the nearby white dwarf LB 9802. 
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Figure fl4.jpg goes here 



Fig. 14. — The circle shows the position of WD0235-12.5 and the neighbouring object, as 
viewed in the IRAC 2 (4.5/im) passband. The separation is approximately 3 pixels (~ 3.6"). 




Fig. 15. — The open points are the flux measurements for WD0235-125 taken with a standard 
3 pixel aperture (to maintain consistency with the bulk of our sample). The solid circles 
indicate the flux measured with a 2 pixel aperture and the filled triangles show the flux 
measured for the companion. The solid line shows a blackbody of 32,400K. The dashed line 
shows a blackbody of 5000 K. 
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Fig. 16. — The left panel shows WD1614+136 (compared to a 22400 K black body) and the 
right-hand panel shows WD1353+409 (compared to a 23600 K black body)/ The dashed 
curve in the right-hand panel shows a 5000 K black body fitted to the point source ~ 5" 
from the white dwarf. Once again, solid points indicate photometry measured in a 3 pixel 
aperture and open circles in a two pixel aperture. 
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Figure fl7.jpg goes here 

Fig. 17. — This shows the zoom in around WD1353+409. The presence of a bright nearby- 
star has caused some banding in this (IRAC1) image, but the star is still well-measured. 
Furthermore, the image reveals a marginal detection of another object ~ 4 ( ~ 5") pixels 
away. 

Table 1. Sample Characteristics 



Name 




logg 


M 


d 








l M ©J 


(pcj 


EUVE0138+25.3/WD0136+251 


39400(200) a 


9.12(13) a 


1.32(8) 


76 


EUVE0237-12.3/WD0235-12.5 


32400(200) a 


8.64(4) a 


1.03(2) 


66 


EUVE0443-03.7/WD0440-038 


65140(200) a 


9.12(12) a 


1.33(7) 


144 


EUVE0521-10.4/WD0518-105 


33000(300) a 


8.70(8) a 


1.07(5) 


99 


EUVE0916-19.7 


56400(2600) a 


9.12(20) a 


1.33(11) 


164 


EUVE1659+44.0/WD1658+441 


30510(200) b 


9.36(7) b 


1.41(4) 


27 


EUVE1024-30.3/WD1022-301 


34800(500) c 


9.20(12) c 


1.27(3) 


61 


REJ0003+43/WD0001+433 


42400(1200) a 


9.30(12) a 


1.37(6) 


101 


REJ0348-005/WD0346-011 


43200(500) a 


9.21(5) a 


1.37(3) 


29 


RE J 1 746-706/WD 1 740-706 


46800(800) c 


8.94(8) c 


1.18(3) 


79 


EUVE0317-85.5/WD0325-857 


40000(10000) 


9.45(15) 


1.33(3) 


36 


WD1614+136 


22400 d 


7.34 d 


0.33 


110 


WD1353+409 


23600 d 


7.54 d 


0.40 


126 


EUVE0653-164/WD0652-563 


35200(600) a 


8.88(10) a 


1.18(6) 


107 



a from Vennes et al 1997 
b from Schmidt et al 1992 
c from Marsh et al 1997 
d from Bergeron et al 1992 



Note. - - The masses are derived from models with hydrogen atmospheres. 
Approximate distances are derived by matching the absolute magnitude derived 
from the spectral fits to the apparent magnitude. 
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Fig. 18. — The solid points show the flux measurements for GD 50. The dotted curve is a 
black body of the appropriate temperature. The dashed curves show the emission from a 
flat, passively reprocessing black body disk around GD 50. The short dashed curve assumes 
an inner disk radius of Oj = 10 10 cm and the long dashed curve has = 10 11 cm. In both 
cases we assume an inclination of 85° to the line of sight and outer radius a Q = lOaj. 
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Fig. 19. — The solid points show the flux measurements for GD 50. The dotted curve is a 
black body of the appropriate temperature. The short dashed curve shows the emission from 
a flat disk composed of superheated dust, as described in the text. Once again, the inclination 
is assumed to be 85°. The inner and outer radii are 10 10 cm and 10 11 cm respectively. The 
long dashed curve shows the level of contribution allowed by the observations if the dust disk 
is optically thin to incident radiation. The inclination in this case is 45° and Tuv — 0.05. 
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Fig. 20. — The dotted line is the white dwarf photosphere emission. The dashed line shows 
the spectrum of a 10 8 years old, 10Mj brown dwarf and the solid line is the combined 
emission. The model is ruled out (at the 2a level) by the photometry at 4.5/im. 
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Fig. 21. — Solid points indicate models of planets from BSL that could be detected by our 
observations if in orbit about WD1614+136. In the case of open circles, these models would 
be undetectable because they lie well below the photospheric emission. The cooling age of 
the white dwarf is only ~ 2 x 10 7 years, but if the planet was formed at the same time as 
the original main sequence star it could be considerably older. 
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Table 2. IR Fluxes of Hot White Dwarfs 



Name 


my 




mj 


F3.6 


F4.5 


F5.8 


F 8 












(mJy) 


(mJy) 


(mJy) 


(mJy) 


H/U v H/U looTZd.o/ vvuuiou-rzoi 


15 


87(3) 


1 P, 


^0 yu ) 


U.U^ 1 1 Zi J 


u.uz 1 yo j 


■f n m 9 


■f n ni 9 


EUVE0237-12.3/WD0235-12.5 


15 


4(1) 


15 


74(6) b 


0.117(2) 
0.107(1) 


0.075(3) 
0.065(2) 


0.043(11) 
0.040(9) 


0.047(10) 
0.031(8) 


EUVE0443-03.7/WD0440-038 


16 


7(2) 


17 


38(5) b 


0.019(2) 


0.022(2) 


< 0.011 


< 0.011 


EUVE0521-10.4/WD0518-105 


15 


82(3) 


16 


76(20) b 
94(2) b 


0.052(2) 


0.037(3) 


0.011(10) 


< 0.011 


EUVE0916-19.7 


17 


3 


13 


0.020(2) 


0.011(3) 


0.018(11) 


< 0.011 


EUVE1659+44.0/WD1658+441 


14 


62 


15 


47(10) b 


0.129(2) 


0.078(2) 


0.039(10) 


< 0.009 


EUVE1024-30.3/WD1022-301 


16 


7(3) 


16 


97(8) 


0.037(2) 


0.024(3) 


< 0.011 


< 0.010 


REJ0003+43/WD0001+433 


16 


8(3) 


17 


0(1) 


0.020(2) 


0.014(3) 


< 0.011 


0.020(10) 


RE J0348-005 /WD0346-0 1 1 


14 


04(2) 


14 


97(1) 


0.238(2) 


0.155(3) 


0.113(12) 


0.044(12) 


REJ1746-706/WD1740-706 a 


16 


6(3) 


17 


2(1) 


0.031(1) 


0.015(1) 


0.009(6) 


0.010(6) 


EUVE0317-85.5/WD0325-857 


14 


8 


15 


43(2) 


0.156(2) 


0.094(2) 


0.050(10) 


0.029(10) 


WD1614+136 


15 


24 


15 


84(2) 


0.116(2) 


0.065(2) 


0.053(10) 


0.014(10) 


WD1353+409 


15 


55 


16 


20(3) 


0.080(2) 


0.052(3) 


0.049(9) 


0.035(9) 


EUVE0653-164/WD0652-563 


16 


4 


17 


4(1) 


0.019(2) 


0.018(2) 


0.023(11) 


0.022(10) 



a The exposure length for this object was twice as long as for the others. 
b taken from Green, Ali & Napiwotzki 2000. 



Note. - - The quoted IRAC errors are statistical. An additional calibration uncertainty of ~ 2% 
(Fazio et al 2004) is applicable. 
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Table 3. Constraints on Disks 



Name 


D 


F 8 


cosi 


T~UV 


M 20 




(PC) 


(mJy) 






(g) 


EUVE0138+25.3/WD0136+251 


76 


< 0.012 


< 0.024 


< 0.10 


< 0.03 


EUVE0237-12.3/WD0235-12.5 


66 


< 0.012 


< 0.07 


< 0.30 


< 0.10 


EUVE0443-03.7/WD0440-038 


144 


0.047 


< 0.03 


< 0.34 


< 0.11 


EUVE0521-10.4/WD0518-105 


99 


< 0.011 


< 0.04 


< 0.16 


< 0.05 


EUVE0916-19.7 


164 


< 0.011 


< 0.10 


< 0.44 


< 0.14 


EUVE1659+44.0/WD1658+441 


27 


< 0.011 


< 0.002 


< 0.010 


< 0.003 


EUVE1024-30.3/WD1022-301 


61 


< 0.009 


< 0.017 


< 0.06 


< 0.02 


REJ0003+43/WD0001+433 


101 


0.020 


< 0.07 


< 0.30 


< 0.09 


REJ0348-005/WD0346-01 1 


29 


0.044 


< 0.013 


< 0.05 


< 0.02 


REJ1746-706/WD1740-706 


79 


0.010 


< 0.022 


< 0.09 


< 0.03 


EUVE0317-85.5/WD0325-857 


36 


0.029 


< 0.013 


< 0.06 


< 0.02 


WD1614+136 


110 


0.014 


< 0.060 


< 0.25 


< 0.08 


WD1353+409 


126 


0.035 


< 0.200 


< 0.843 


< 0.265 


EUVE0653-164/WD0652-563 


107 


0.022 


< 0.09 


< 0.37 


< 0.12 



Note. — Constraints on inclination for a flat, blackbody disk and on mass for 
an optically thin dust disk. 
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